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Particle Correlations at RHIC - Scrutiny of a Puzzle 
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We present calculations of two-pion and two-kaon correlation functions in relativistic 
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heavy ion collisions from a relativistic transport model that includes explicitly a first-order 
phase transition from a thermalized quark-gluon plasma to a hadron gas. We compare 
the obtained correlation radii with recent data from RHIC. The predicted -R S ide radii agree 
with data while the R out and -Riong radii are overestimated. We also address the impact 
of in-medium modifications, for example, a broadening of the p-meson, on the correlation 
radii. In particular, the longitudinal correlation radius -Ri on g is reduced, improving the 



comparison to data. 



Two-particle correlations at small relative momenta have been predicted to be partic- 
ularly sensitive to a phase transition from quark-gluon matter to hadronic matter [|TJ. 
For a first-order phase transition, larger hadronization times were expected to lead to 
considerably enhanced correlation radii, characterizing the space-time extension of the 
particle-emitting source, compared to, for example, a purely hadronic scenario. The radii 
should also depend on the critical temperature T c , the latent heat, the initial specific 
entropy density or the initial thermalization time of the quark-gluon phase. 

Here, we discuss relativistic transport calculations at RHIC energies that describe the 
initial dense stage by hydrodynamics M and the later more dilute stages by microscopic 
transport || of the particles. The two models are matched at the hadronization hyper- 
surface [§. In the hadronic phase the particles are allowed to rescatter and to excite 
resonances based on cross sections as measured in vacuum. For the initial dense (hydro- 
dynamical) phase a bag model equation of state exhibiting a first-order phase transition 
is employed. Hence, a phase transition in local equilibrium that proceeds through the 
formation of a mixed phase, is considered. The details of this relativistic hybrid transport 
model can be found elsewhere ||. 

We first briefly summarize the main conclusions obtained in previous work |5],|6],|7|] . Then, 
we show in detail the results for pions when calculating explicitly the correlation functions 
from the source function of the transport model. These results are subsequently compared 
to calculations that take in-medium modifications into account and to experimental data 
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for central Au+Au collisions at y/s NN = 130 GeV |||| . 

Model calculations have demonstrated that the freeze-out hypersurfaces of pions extend 
to rather large radii and times compared to the size of the mixed phase []7|||. In this 
hadronic phase many soft collisions take place which hardly modify the single-particle 
spectra but have a strong impact on the correlation functions that measure the final 
freeze-out state. Studying this in detail leads to the following conclusions J7J: (i) The 
dissipative hadronic phase leads to a rather large duration of emission, (ii) The R on t/Rside 
ratio, thought to be a characteristic measure of this emission duration, increases with 
transverse momentum, (iii) The specific dependencies of the interferometry radii on the 
QGP properties are rather weak due to the dominance of the hadronic phase. This even 
results in qualitative differences if calculations with and without this subsequent hadronic 
phase are compared (dependence on the critical temperature) 0. 

The correlations of kaons provide a severe test of the pion data and have several advan- 
tages 0. In particular, the kaon density is much lower than the pion density |1(J. Hence, 



multiparticle correlations that might play a role for the pions are of minor importance 
for the kaons. Also, the contributions from long-lived resonances are under better control 
for kaons. The i? ut/-Rside ratio for kaons is shown in Fig. 1. The sensitivity to T c and 
to the specific entropy density (SPS vs. RHIC) is enhanced at larger transverse momenta 
(Kt ~ 1 GeV/c). This enhanced sensitivity is also driven by a strong increase of the di- 
rect emission component (from the phase boundary) at high Kt as shown in Fig. 2. More 
and more kaons (up to ~ 30%) escape the initial stages (unperturbed by the hadronic 
phase). 
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Figure 1. R out / Rside for kaons at RHIC 
(full symbols) and at SPS (open symbols), 
function of Kt for critical tempera- 
tures T c ~ 160 MeV and T c ~ 200 MeV, 
respectively. 



Figure 2. Fraction of kaons r or i gin that ori- 
gin from a particular reaction channel prior 
to freeze-out. These are resonance decays 
(full circles), direct emission from the phase 
boundary (stars), elastic meson- meson (dia- 
monds), or elastic meson-baryon (open cir- 
cles) collisions. 
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Figure 3. Pion correlation parameters R a , R s , R\, and A as obtained from the 
'QQP+hadronic rescattering' transport model calculations + correlation after burner (T) 
compared to experimental data from PHENIX and STAR. For theory (T), the effects 
of a finite momentum resolution (f.m.r.) and a broadened in-medium p spectral function 
(pi.m.) are calculated. Note the strong reduction of Ri when an in-medium p is considered. 

The pion correlation parameters R Q , _R S , i?iand A are obtained from the explicit cal- 



culation of the correlation functions |L1[ in the respective transverse momentum bins 
and subsequent fitting of these correlation functions to a Gaussian form of the correlator 
C2 = 1 + Aexp(— R^ql — RgQs ~ R\Qi)- The results of fitting the 3-dimensional pion cor- 
relation functions (as obtained from the transport calculations (RHIC initial conditions, 
T c « 160 Me V) + correlation after burner (by Pratt) JlIJ) are shown, together with the 
experimental STAR and PHENIX data in Figure 3. Even for a first-order phase 

transition scenario, the calculated interferometry radii (T) are not unusually large. Fi- 
nite momentum resolution (f.m.r.) reduces the true correlation radii and the A intercept 
parameter, in particular, at higher 0. This is corrected for in the experimental anal- 
ysis. While the -R S id e radii show agreement the R out radii are too large compared to the 
data. However, these pion radii are considerably smaller than the corresponding radii 
obtained from the coordinate-space points and using expressions for the Gaussian radius 
parameters based on a saddle-point integration over the source function p||6|,|T|,|r2| . Only 
the calculations as shown in Figure 3 which are obtained from the complete calculation 
and the performed fits should be compared to data. The i? out /-Rside ratio is also larger 
than unity for the fitted values and confirms the so-called HBT-puzzle, i.e., the RHIC 
data from STAR and PHENIX HU indicate a decreasing ratio with K T (even below 1) 
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and the calculations show only ratios larger 1. We remark that at the CERN-SPS (see 
Refs. in |J) the -R ou t/-Rside ratio seems to increase with K T , being larger than 1, im- 
plying a real qualitative change of the reaction dynamics from SPS to RHIC energies. 
A discussion of theoretical assumptions and uncertainties, experimental corrections, and 
possible solutions of this problem (including strongly opaque sources or scenarios with 
strong supercooling leading to spinodal instabilities) is provided in 0. 

In-medium modifications may alter the correlation radii due to their impact on the 
mean free path via modified cross sections at finite temparatures or densities (Tf. G ~ 
100 — 120 MeV). Increasing the width of the p from the vacuum value of T p = 150 MeV 
to T^ m ' = 350 MeV enhances the opacity of the system (more tttt «-> p processes), with a 
particular reduction of the radius -Ri ong , improving the comparison to data (Fig. 3). 

We have demonstrated that the kaon interferometry measurements, in particular at 
high Kt, will provide an exellent probe of the space-time dynamics (close to the phase 
boundary). In addition they represent a severe test of the pion correlations and may help 
to better understand the HBT-puzzle flI3| , i.e., the difference in the Kt dependence of the 
i? out /-Rside ratio between the model predictions and the experimental RHIC data. The 
differences are due to the R ont radii (which are larger in the model calculations) while the 
.Rside radii seem to be described reasonably. An in-medium broadening of the p-meson 
reduces in particular -Ri on g, improving the comparison to data. 
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